Nematic and Polar order in Active Filament Solutions 
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Using a microscopic model of interacting polar biofilaments and motor proteins, we characterize 
the phase diagram of both homogeneous and inhomogeneous states in terms of experimental pa- 
rameters. The polarity of motor clusters is key in determining the organization of the filaments in 
homogeneous isotropic, polarized and nematic states, while motor-induced bundling yields spatially 
inhomogeneous structures. 
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Soft active systems are exciting examples of a new 
type of condensed matter where stored energy is contin- 
uously transformed into mechanical work at microscopic 
length scales. A realization of this are polar filaments 
interacting with associated molecular motors in the cell 
cytoskeletoni*^. These systems arc characterized by a 
variety of dynamic and stationary states which the cell 
accesses as part of its cycle2i4i^. 

There have been a number of recent theoretical stud- 
ies of the collective dynamics of isotropic and polar- 
ized solutions of active filaments. These include numer- 
ical simulations^, mesoscopic mean-field kinetic equa- 
|-jpj^g7j9ji0ji^^ and hydrodynamic equations where the sys- 
tem is described in terms of a few coarse-grained fields 
whose dynamics is inferred from symmetry considera- 
tionsiSiiiiiii^. Previous work has focused on how mo- 
tor activity renders homogeneous states unstable to the 
formation of spatial structures, such as bundles, vortices 
or asters. In this article we study the profound effect of 
motor activity on the possible homogeneous states of the 
systemic. We also find important differences in the na- 
ture of the instabilities from these homogeneous phases. 
Starting from a microscopic model of interacting rigid fil- 
aments, we characterize a phase diagram of homogeneous 
and inhomogeneous states of active filaments in terms of 
experimentally variable parameters. We find in partic- 
ular that the formation of a non-equilibrium polarized 
phase at low densities can be driven by motor polarity 
without the need for filament polymerization^Sii^. 

We describe the system by a concentration of polar 
filaments /(r, n,t) in two dimensions [d = 2), modeled 
as hard rods of fixed length £ and diameter b {£ >> b) 
at position r with filament polarity characterized by a 
unit vector n, and a density of motor clusters m(r, t) at 
position r. The filament and motor concentrations satisfy 
the equations 



ume, J^{r,h,t) = --f^f{r,h,t)'RVx{r,n,t) , with 



dtf 
dtm 



-V • J™ 
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where TZ ~ h x dn and the translational (J/,J„i) and 
rotational (JT) currents have diffusive, excluded volume 
and active contributions. The rotational current is JT" = 
+ + J^, with a diffusive current J^{r,n,t) = 
—DrTZf{r,h,t) and a contribution from excluded vol- 



Vx{r,ni,t)^kBT 



|ni X n2|/(r + |,n2,t) 



(3) 

where ^ = fiisi — n2S2 is the separation of the centers 
of mass of the two filaments and ... = J^£/2 de- 
notes an integration along the length of the filament, 
parametrized by s. The active contribution to the ro- 
tational current (low density approximation) is 

J Si J S2 ^ h2 

x/(r,ni,i)/(r + |,n2,i) , (4) 
where the motor induced angular velocity is written as 

o;-^ = 2(70 + 7ini • 62) (ni X n2) . (5) 



It consists of two parts, corresponding to two classes of 
motor clusters (see Fig. : polar clusters which tend to 
bind to filaments with similar polarity ( 70/71 3> 
and non-polar clusters which bind to filament pairs of 
any orientation ( 70/71 <IC 1)^. Earlier work by two of 
us^ considered only non-polar clusters (70 = 0)ii. Both 
7o and 71 will increase with increasing binding rate of the 
clusters to the filament. A passive polar crosslink (e.g. 
a-actinin on F-actin) will also have this effect. Since 
the binding rate can have different dependence on ATP- 
consumption than the motor stepping speed, we expect 
a different dependence on ATP hydrolysis rate than the 
active contributions to the translational currents. 

The translational currents are J / = + + J-^ and 



-AnVm + x / / n/(r, n, i)m(r + ns,t) 

'J s 'J h 



(6) 



where x depends on the speed and processivity of the 
motors. The filament diffusive current, JP = —Dijdjf, 
is expressed in terms of a diffusion tensor Dij = (Z?|| + 

Dj_)S,,/2 + iD\\-Di_)Q,, withQ - - - ■ ^ 



2 5ij . The ex- 



riiTij- 

cluded volume contribution is J/^ = —-^j^fdjVx- The 
active contribution to the translational current is 

J Si J S2 ^ n.2 

x/(r,ni,i)/(r + ^,n2,i) , (7) 
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FIG. 1: Polar and nonpolar clusters interacting with polar 
filaments. Assuming that clusters always bind to the smallest 
angle, polar clusters {g — > oo) bind only to filaments in con- 
figuration (a) while non-polar clusters {g = 0) bind to both 
configurations equally. 



with V the motor-induced velocity, taken of the form 



v"^ = |(n2 - ni) + - A(ni + fia) 



a ^ 
22i 



(8) 



The parameters a, (3 and A have dimensions of veloc- 
ity and depend on the angle between the filaments. The 
term proportional to (3 drives the separation of filaments 
of opposite polarity, while the A contribution arises from 
the net velocity of the filament pair—. The negative sign 
reflects the fact that filament mean motion due to motor 
activity is opposite to their polarity. The contribution 
proportional to a arises from spatial variations in motor 
activity along the filament, such as motors stalling before 
detaching at the polar end. It drives bundling of fila- 
ments of the same polarity. These parameters where esti- 
mated in Ref— via a microscopic model of motor- induced 
filament dynamics as /3 ~ A ~ uq, a ~ uo{lm/l) << uq, 
with uo the mean motor stepping rate and Im the length 
scale (of order of the motor cluster size) for spatial varia- 
tions in motor activity. As seen below, this term is crucial 
for developing inhomogeneities and pattern formation. 

To study the macroscopic properties of the solution, 
we truncate the exact moment expansion of /(r, fi, t) as 



mean value (we let itiq = mb"^) and the filament den- 
sity, /(r, n,t), and its moments can be expressed in 
terms of their spatial averages, i.e. -j J drSij{r,t) = 
Sij{t), ^ J (irp(r,t) = p{t), with A the area of the sys- 
tem. In the following all lengths are measured in units 
of the filament length, i. Averaging over the orientation 
n using Eq. , we find that in a homogeneous system 
P — Po = constant and 



^^^"^ mo(27o-7i) paSijPj 



dtSij 



3tt 

4:Dr - 



(11) 
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+2mQpolo{piP] -\5ijp^^ . (12) 



In a passive system (70 = 71 = 0) there is a transition 
from an isotropic state to a nematic state. A mean-field 
description of such a transition, which is continuous in 
2d (but first order in 3d), requires that one incorporate 
cubic terms in the nematic order parameter in the equa- 
tion of motion. The transition here is identified with the 
change in sign of the decay rate of 5"^, which signals 
an instability of the isotropic homogeneous state. This 
occurs when excluded volume effects dominate at a den- 
sity pn — 37r/2. The homogeneous state is isotropic for 
Po < Pn and nematic for po > Pn- No homogeneous po- 
larized state with a nonzero mean value of p is obtained 
in a passive solution. 

We now turn to an active system. We introduce a di- 
mensionless filament density, p ~ Po/ Pn, a dimensionless 
motor cluster activity, p = pisirriQ^i/ D^, and a parame- 
ter measuring the polarity of motor clusters, g = 70/71 
with .9 = corresponding to non-polar clusters. Time is 
measured in units of . The steady states of Eqs. (|11|) 
and H12|l are the stable solutions of 

= -aipi + bipSijPj , (13) 

= -a2Sij + b2p(piPj - , (14) 



/(r, n, t) = { 1 + 2p(r, t)-h + 45,, (r, t)Q,, } , (9) where 

keeping only the first three moments, 

J dh /(r, n, t) = p{r, t) (density), 

j dh h f{r,h,t) ~ p{r,t)p{r,t) (polarization), (10) 

J dh Qij /(r, n, t) = p(r, t)Sij (r, t) (nematic order) . 

a. Homogeneous Bulk Steady States. We first con- 
sider the dynamical equations for a spatially homoge- 
neous solution. In this case the only contributions to 
the equation of motion of the filament density come from 
rotational currents. The motor density has a constant 



fli = 1 - pgp , 

a2 = 4[l-p(l + /i/4)] , 



(15) 
(16) 



and bi = A+p (1 — 2g), &2 = 2(7 /i. At low density the only 
solution is Pi = and Sij ~ and the system is isotropic 
(I). The homogeneous isotropic state can become unsta- 
ble in two ways. As in the passive case, a change in 
sign of the coefficient 02 signals the transition to a ne- 
matic (N) state. Motor activity lowers the density for 
the I-N transition which occurs at pin{p) = 1/(1 + A'/4). 
At p > pin{p) the solution acquires nematic order, 
with = So{ninj — Sij/2), where the unit vector n 
denotes the direction of broken symmetry. We obtain 
an expression for the amount of nematic order (Sq) by 
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adding a cubic term —C2p^ SuSkiSij to Eq. ifT^ giving 
Sq = i-^— 202/02. The isotropic state can also become 
linearly unstable via the growth of polarization fluctua- 
tions in any arbitrary direction. This occurs above a sec- 
ond critical filament density, pip{^) = l/{gfj), defined 
by the change in sign of the coefficient oi controlling the 
decay of polarization fluctuations. For p > pip(fj,) the 
homogeneous state is polarized (P), with pi ^ 0. The 
alignment tensor also have a nonzero mean value in the 
polarized state as it is slaved to the polarization. One 
can identify two scenarios depending on the value of g. 

I) For g < 1/4, the density pjp is always larger than 
Pin and a region of nematic phase exists for all values 
of jj,. At sufficiently high filament and motor densities, 
the nematic state also becomes unstable. Fluctuations in 
the alignment tensor are uniformly stable for 02 < 0, but 
polarization fiuctuations along the direction of broken 
symmetry become unstable for oi < pbiSo/2, i.e., above 
a critical density 




FIG. 2: The homogeneous phase diagram for g < 1/4 (the 
figure is for g — 1/10 and C2 = 50). The vertical axis is po/pjv, 
the horizontal axis is p. 
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2c2i?(l - R) 



C2R 



(17) 



where R — pin / Pip- 
has = Vrini and 



The polarized state at p > pNP 



Pi = Pqu^ 

_ 2aia2_ fi 

i2 L-"- 



^-^i and 



FfcTfeL- yifkYl ana = ^0^1 " 'l^P^ = 
2|ai/6i|. The "phase diagram" is shown in Fig.|21 

II) When g > 1/4, the boundaries for the I-N and 
the N-P transitions cross at px = 1/(3 — 1/4) where 
Pin = Pip = Pnp and the phase diagram has the topol- 
ogy shown in Fig. 13 For p > p^ the system goes di- 
rectly from the I to the P state at pip, without an in- 
tervening N state. At the onset of the polarized state 
the alignment tensor is again slaved to the polarization 
field, poSij = ^ {prPj - \5ijp^) , and pi = poUi and po 
is determined by cubic terms in Eq. ITTI 

If 71 = 0, with 7o 7^ 0, the I-N transition is indepen- 
dent of motor density and always occurs at po = pN- 

b. Spatially inhomogeneous states. Spatial inhomo- 
geneities of course affect the stability of the homogeneous 
states described above. As shown by several authors, 
the rate of motor-induced filament bundling can exceed 
that of filament diffusion yielding the unstable growth 
of density inhomogeneities^iSiiSiii. States with spatially 
varying orientational order, where the filaments sponta- 
neously arrange in vortex and aster structures, are also 
possibl9iSiiS4i&. To examine the role of spatial inhomo- 
geneities, we have obtained coupled equations for the first 
three moments of the filament concentration defined in 
Eq. © by an expansion in spatial gradients described 
elsewhereii*iS. These equations can then be used to 
study the linear stability of each of the homogeneous 
states against the growth of spatially- varying fluctuations 
in the hydrodynamic fields. These are the fields whose 
characteristic decay time are much longer than any mi- 
croscopic relaxation time and become infinitely long lived 
at long wavelengths. We find that the low frequency hy- 



Sp{ninj 
So 



5ij/2), where 



1 



Inhomogeneous 




FIG. 3: The phase diagram for g > 1/4. For p > px, where 
Pin and pip intersect, no N state exists and the system goes 
directly from the I to the P state. Inhomogeneous states form 
for po > pb- The ps line may lie above the pjvp — pip line, 
as shown in the main figure {g = 1, 71/a — 2.3), or cross 
through the N and I states, as shown in the inset {g = 1, 
71 /a = 0.6), depending on the value of ji/a. (with C2 — 50). 



drodynamic modes of this active system are determined 
by fluctuations in the conserved densities and in variables 
associated with broken symmetries. A change in sign in 
the decay rate of these modes signals an instability of 
the macroscopic state of interest. For simplicity we only 
discuss here the case of constant motor density. 

Isotropic state. This has been studied in Re^ii. The 
only hydrodynamic variable is the filament density, p. 
The decay rate of Fourier components oi Sp ~ p ~ po 
at wavevector k is controlled by the interplay of diffu- 
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sion and motor-induced filament bundling described by 
a. The homogeneous I state is unstable at large length 
scales for po > pb, with pB ^ D\\/(mQa) ~ ji/{ap). 
The homogeneous state is stabilized at short scales by 
excluded volume corrections and higher order terms in 
the density gradients. The density instability is driven 
entirely by the bundling rate a, while /3 plays no role. 

Polarized state. The hydrodynamic modes in the P 
state are those associated with fluctuations in the fila- 
ment density and in the director field, n{r,t), defined by 
p{r,t) = p(r, t)n(r, t), with |n| = 1. The coupled hy- 
drodynamic modes describing the decay of Fourier com- 
ponents of density, 6p and director fluctuations, dn — 
n — y = xSux of wavevector k are always propagating, 
with velocity whose magnitude and direction are con- 
trolled by both the activity parameters /3 and A. For k 
along the broken symmetry direction, the modes decou- 
ple (i.e., 6p ~ e^p*-*^^*, i5na; e^"'-'^^*) and are given by 



fc2 

ikcipp(3 — 



1 — ^ — 20ppa 



(18) 



5^2 2 -1 
z„ = -ikc2pp(3 ~ ^ [l + -^PKg - 6a)J , (19) 

where f3 = (3/ji, a = a/71, and ci and C2 are numbers of 
order one. We have used = and A ~ /3. Like 

the I state, the homogeneous P state is linearly unsta- 
ble for pq > Pb- The nature of the instability changes, 
however, from diffusive in the I state to oscillatory in 
the P state, suggesting that spatially inhomogeneous os- 
cillatory structures, such as vortices, may be stable at 
high fllament or motor densities. The rotational effects 
described by ^ ^ 71 stabilize director fluctuations, but 
destabilize the density. 

Nematic state. The hydrodynamic variables in the 
N state are again the fllament density and a director 
field n(r,t), defined in terms of the alignment tensor as 



SoiniUj 



The decay of density and direc- 



modes which are always diffusive. The modes decouple 
for k along the broken symmetry direction, with 



Zp — k 



i - 2ppa 

19 

+ 36 



pp 



(20) 
(21) 



tor fiuctuations is controlled by coupled hydrodynamic 



Once again, the homogeneous N state is destabilized by 
the growth of density fluctuations for pQ > ps, while long 
wavelength director fluctuations remain stable. For ar- 
bitrary direction of k relative to the direction of broken 
symmetry director fluctuations also become unstable at 
high density, but the fastest growing mode is always as- 
sociated with the build-up of density inhomogeneitics. In 
contrast to the P state, the instability is always diffusive. 

To summarize, we have studied the phase behavior 
of polar fllaments interacting with polar clusters. We 
have shown that in addition to the homogeneous isotropic 
phase, both homogeneous polarized and nematic states 
can be obtained as a function of filament density and mo- 
tor activity and polarity. The instabilities of the homo- 
geneous states are controlled by the bundling rate a and 
occur for po > pB ^ 7i/(a/-i). The location of this line 
in the phase diagram depends crucially on the parameter 
7i/a. If 7i/q! > l/g ~ 71/70, then pB > pp as shown in 
Fig. 13 and the homogeneous nematic state is always sta- 
ble, when it exists. 1171/0; < l/g = 71/70, then pB < Pp 
as shown in the inset of Fig. and all homogeneous states 
can be destabilized by fllament bundling, albeit through 
different (diffusive versus oscillatory) mechanisms. 
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